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In order to estimate the effect of dynamical gluons to chiral condensate, the gap equation of 5'L'^(2) 
gauged Nambu-Jona-Lasinio model, under a constant background magnetic field, is investigated up 
to the two-loop order in 2 + 1 and 3 + 1 dimensions. We set up a general formulation allowing 
both cases of electric as well as magnetic background field. We rely on the proper time method 
to maintain gauge invariance. In 3 + 1 dimensions chiral symmetry breaking (xSB) is enhanced 
I by gluons even in zero background magnetic field and becomes much striking as the background 

. field grows larger. In 2 + 1 dimensions gluons also enhance x^B but whose dependence on the 

background field is not simple: dynamical mass is not a monotone function of background field for 
a fixed four-fermi coupling. 
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, The Nambu-Jona-Lasinio (NJL) model the model of four fcrmion interaction, has been discussed as a possible 
' realistic mechanism of chiral symmetry breaking (xSB) |^,^ . Although the interaction is nonrenormalizable in more 
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I. INTRODUCTION 



than l-fl spacetime dimensions, it is regarded as a low energy effective theory of elementary fermions after integrations 
of Higgs and/or gauge fields, and has been used as a phenomenological model to describe hadronic spectra, decays, 
T-H _ and scattering (While these properties of hadrons should be derived from quantum chromodynamics, analytical 
. methods on hand are so limited that dynamics in the low-energy region could not be easily explored.) 

In order to grasp qualitative behavior of a system, it is sometimes useful to examine response to external fields 
or sources. Many such attempts have been made so far, for example, the 0(3) gauge-Higgs model in a magnetic 
field source fermionic models minimally coupled to strong electromagnetic fields Q, the NJL model minimally 
coupled to a constant magnetic field and curved spacetime 0] and its extension to supersymmetric NJL model 
i-S^ the instanton motivated four-point interaction model of fermion at finite density and so forth. 

I Among these, the most interesting outcome is found in the study of the NJL model minimally coupled to constant 
. external magnetic fields pO|-p^; where mass generation occurs even at the weakest attractive interaction between 
/\ ' fermions in terms of so called the "dimensional reduction(DR)" [Q. The case is also investigated by the present 
authors to find that the origin of DR is the infrared divergences followed from the fermion loop integral under the 
influence of the external magnetic fields . The phenomenon is now understood universal provided interactions are 
shortrange |Q. 

Inclusion of dynamical gauge fields is made for QED, in terms of the Schwinger-Dyson equation [ p^ or renormal- 
ization group (RG) [p^ . The results show that the dynamical symmetry breaking always occur with the aid of an 
external magnetic field. 

The motive for this work is traced to that of Gusynin, Hong, and Shovkovy where 2 -I- 1-dimensional SU{2) 
gauge theory is investigated in a constant magnetic field by using a constant gauge potential with translational 
invariance |19|| , to find that magnetic catalysis of x^B does not occur. The study in 3-1-1 dimensions is made in the 
reference ]18[|, where SU (N) gauged NJL model is handled in the case of weak as well as constant magnetic field: 
they incorporate dynamical effect of gluons in RG to reach the conclusion that existence of the external field does not 
change the condition of xSB. Each result indicates, contrary to our expectation that gluons always trigger xSB (since 
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they do in the low energy region), that dynamical gluons in external magnetic field do not play a major role to xSB. 
This curious situation is our starting point. There needs a more detailed analysis. 

In this article we study xSB under a constant background magnetic field in SU{2) gauged NJL model, paying 
attention to (i) the direct effect of dynamical gluons to the gap equation, not in terms of RG |l7yi^ , and (ii) the 
effect of gauge choice to the results: since in nonabelian case the situation, = B (constant > 0) ; others = 0, the 
choice [l^] 

AI^AI^VB; others = , (1) 

cannot connect with [ p^ , 

B B 

Al^--y, Al = -x, others = : (2) 

since the remnant gauge transformation which leaves B^ invariant, that is, any gauge transformations with respect 
to the third axis cannot bring the gauge (|l|) to the gauge (|^) or vise versa. To work with the Green's function in 
momentum space it is convenient to adopt eq.(0), but this is not a solution of the nonabelian Maxwell equation in 
the vacuum, 

[D^F^^r ^ d,F;^ + e^^^^F^^, = . (3) 

In what follows we rely upon the WKB semi-classical approximation where the classical solution plays a fundamental 
role. Therefore we should work with the choice eq.(||), satisfying eq.(^). Moreover the choice (||) fulfills a covariantly 
constant condition, 

{D,F,,T ^ dpF;^ + e^'^M^F^^, = , (4) 

enabling us to use the Fock-Schwingcr proper time method that was originally developed to handle abelian 
background fields and that gives us a gauge independent result. 

For calculations of the functional determinant, a regularization for ultraviolet divergence is necessary. We introduce 
a proper time cutoff, which preserves gauge symmetry. In addition, there needs another cutoff: in order to estimate 
the contribution from dynamical gluons, two-loop order of the effective potential must be taken into account. The 
gluon propagator under the magnetic background suffers from the famous instability, due to tachyonic modes of the 
spin 1 propagator under the background magnetic fields [ pT[ . To circumvent this instability we introduce a gluon 
mass whose square is larger than the magnitude of a magnetic field. 

The paper is organized as follows: in Sec. || a general formulation is present ed, where we assume generic backgrounds 



satisfying the covariantly constant condition eq.(M). In the next sections. Sec. Ill and Sec. IV, the gap equations in 2-1-1 
and 3 -f 1 dimensions are shown. The last section Sec. ^ is devoted to discussion. In Appendix ^ calculations of the 
kernel is presented and in Appendix ^ the gluon propagator is represented by means of the proper time method. Then 
in Appendix ^ we make an explicit proof that our classical solution does satisfy the covariantly constant condition. 
And finally in Appendix]^, in order to ensure the gauge independence of our calculations, that is, the correctness of 
those, we study the Ward-Takahashi relation of vacuum polarization function. 



II. FORMULATION 



In this section, we derive the effective potential of SU{2) gauged NJL model. The Lagrangian of the system in the 
Euclidean formulation is given as 



D = 4, 
£» = 3 , 



(5) 



where Af^ = A'^Ta with T^'s (a — 1,2,3) being the SU{2) generators given by Ta = cra/2, where Ca's are the Pauh 
matrices. For the 2 -I- 1-dimensional case, a spinorial representation of the Lorentz group is given by two-component 
spinors, so that corresponding gamma matrices are 2x2. There is no chiral symmetry. In order to be able to discuss 
chiral symmetry, we introduce an addtional flavor such that 



[111) : ^ ^^3 ^ ( V^i, -V^2 ) ^ ( ^1^3, -V4'^3 ) , 



(6) 
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with the 4x4 gamma matrices 



7m 



_ I 
-a. 



; ^ = 1 ~ 3 , 74 



1 

1 



75 = 71727374 = 



il 
il 



Chiral symmetry is realized as 



(7) 



(8) 



yielding a global U{2) symmetry which is broken by a mass term into U{1) x U{1). 
The partition function of the model is read as 



D [gauge] D [ip] D [tp] exp 



d^xC 



D [gauge] D[cr]D[n]D[^/j]D[ip] 



X exp 



where the auxiliary fields, a and tt, have been introduced to erase the four-fermi interactions, 



TT r 



7r75 for I? = 4 . 

7ri74 + 7r275 for I? = 3 



(9) 



(10) 



and a measure of gauge fields, I?[gauge], is specified after the following procedures: (i) first, integrate with respect to 
fermions to give 



Z = D [gauge] D[a]D[K] exp 



LnDet 



l^,{^^,-iA^) + {a + in • F) 



(11) 



(Here and hereafter Ln, Det, and Tr designate the functional logarithm, determinant, and trace, respectively.) (ii) Sec- 
ond, set up an ansatz; a{x) = m{: constant), tt = 0, under which the equation of gauge fields reads. 



If we take Fj^^^ — F^^ = ; F^^, — constant , then 



A^^^O; (1 = 1,2) 



43 f3 



(12) 



(13) 



the equation ( |12| ) is fulfilled. (It is easy to convince that the left-hand-side vanishes but the proof of the right-hand- 
side is rather lengthy, so it is relegated to Appendix ^.) Therefore our classical solution (p^ obeys the nonabelian 
Maxwell equation (p^). We do not take the higher powers of Sa and Stt, 



<t{x) = m + da{x) , 7r(a;) = + S7r{x) ; 



(14) 



but remain the lowest part. In this section we do not restrict ourselves in the pure magnetic case but in generic cases 
where both electric and magnetic backgrounds coexist. (iii) Third, expand the gauge fields around 

= -4^ + ' (15) 
with being designated as quantum fields. Here, with the aid of the Faddeev-Popov trick, the measure is defined 



^[gauge] = D [g;] 



^ SG" 








exp 


262/ 



with 



(16) 



(17) 



The gauge transformation now reads. 
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(18) 



with D^f = 5 df^ + t {A'^^ + Q^), so that the Faddeev-Popov determinant is given by 



(19) 



(In what follows, however, it is not necessary to worry about the F-P terms; since they are irrelevant to the gap 
equation.) 

The partition function is given, up to 0{Q^), by 



Z[A] = exp 

D [Ql] exp 



- / ^^^2;— + LnDet [7^(9^ - aj^Tg) 



X exp 



' 4g2 ■^A'J'-^M"^ 



2e 



where A ^ is the inverse of the gluon propagator. 



^- J d^xd^y {^.-l + e^n^, + e2(F-P term))"^" QI 



(A- 



-l\ah 



2\ab 



1 



(20) 



(21) 



under the gauge (13). In eq.(pTl), the symmetric matrix e'^'^'^T^^, has negative eigenvalues after the diagonalization 
pl| , thus there are tachyonic singularities. This is due to a large magnetic moment of spin 1 particle. (Recall that we 
are in the Euclidean world so that all gauge fields are considered as magnetic one.) In view of eq.(|2^), these tachyonic 
singularities become harmless if we introduce a gluon mass Mg that should obey 



(Mg)2 > 



The term in eq.(|2 
52 



jjafc ^ _ 



-LnDet 



2=0 



TrLn 



(22) 



(23) 



2=0 



is the vacuum polarization tensor. We omit the isospin index 3 of as well as J-'^^ from now on. 
Integrating with respect to the quantum field Q^, we obtain 



Z — exp 



X (m — independent terms) , 



where 



with 



V 



2<?2 



V 



(D) 



V 



(D) 



VT 



-TrLn 



V 



(D) _ 



J^^Tr(nA) = -J^^ J d^xd^y IVf,{x,y)^%{y,x) , 



(24) 

(25) 

(26) 
(27) 



being the one-loop and the two-loop effective potential respectively. Here V is the {D — l)-dimensional volume of the 
system, T is the Euclidean time interval. The stationary condition for the effective potential, 



gives a gap equation. 



(4^) 



D/2 



dm 



= f[^\x) + fi-\x) 







m 



x=j^, (0 < X < 1) 



(28) 



(29) 
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where A is the ultraviolet cut-ofF, and 



J^ \ ) - 2AD-2 Q^2 



i = 1,2 



(30) 



The remaining task is therefore the calculation of the effective potential eqs.(^6|) and (|27|). Let us start with the 
one-loop part. We utilize the proper time method |20|: 



1 



VT 
1 



2VT 



-TrLn 



1 



dT r-^e-"nr / x{x\e-^"°\x) , 



2VT 7i/A^ 

where the ultraviolet cutoff A has been introduced, tr is taken only for the spinor and the isospin indices, and 

1 



i/n = n 



with 



[Xf^,p^] = iS^iy. Write the kernel as 



to find 



X(x,y;r)= (x|e-^^°|2/) 



K{x,y\T) 



(47rr)^/2 



det 



sinrJ^/2 
tT/2 



X exp [- -{x-y)^ 



— cot 

2 2 



'1/2 



(Xo(r)I + if3(T)T3) exp[iT3C] 



whose derivation is given in Appendix In eg. (^) , quantities are defined such that 



1 

2^ 



it'o (t) = cosh — ^ cosh — ^ 75 sinh — ^ sinh — ^ , 



Ksir) = cr^^ f A^+^ sinh — - cosh — - + N^^ cosh - 



where 



F+ = \fW+W 
F_ = 

{\B + E\ + \B - E\} 



R 



B = Ti2]E= {Ti3,T2:i) ; for = 3 
E= {Tii, T2i, J^zi) \ 



{\B + E\-\B-E\} 

^- = 7, I -° = {■'^23, -r 31, ^12) , 



for £> = 4 



and 



N. 



N+ = 5^ = N 



TV" = ; for D = 3 



F2 -F2 

with = e^i,\pT\p/2 being the dual of T^^. 



N, 



p p p p 

liv = p2 _ p2 



for £> = 4 



(31) 

(32) 

(33) 
(34) 



(35) 

(36) 
(37) 
(38) 

(39) 
(40) 



(41) 
(42) 
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Therefore 



1 



pOQ n A poo 

/ drr-ie-™\r/d^xX(x,a;;T) = — — ^ / dr T-^'/^-'e-^"'" GoirF) 
Ji/A^ J (47r)^/^ 7i/A2 



(43) 



where 



■ coth 



GnirF) = 



tF. 



T^F+F^ 



coth 



tF. 



coth I 



for £1 = 3 



for D = 4 



(44) 



In order to calculate the two-loop contribution eq.(g7|) first we express the gluon propagator eq.(|2l|), in terms of 
the proper time, as 



dre 



{x\e 



\y) 



a,6 = 1,2,3 



(45) 



where Mg is the gluon mass eq.(p2|). (Gluons would be massive in the considering situation, that is, in the confining 
phase.) The results whose derivation is relegated to Appendix H are 



^%{x,y) = [(cosC + esinC) {/^l,{x - y) + e /\%[x - y))]' 



= 1,2 , 



-tM„ 



dr 



(47rr)^/2 



exp 



4r 



i^-yf 



(46) 
(47) 



and others — 0. (If we work with an abelian gauge theory, only A^^ term survives. In this sense, we call eq.(p7|) the 
abclian part.) Here C has been given in eq.(Pq), e is a 2 by 2 antisymmetric matrix (e^^ = — e^^ = 1), 



A;,,(a;-2/) = / dr 



(4^r)^/2 
1 



(cos2tJ^)^^ 



det 



sin tJ- 



-, -1/2 



X exp 



[x-y]- {TcoXtT) ■ (x-y) 



(48) 



A^.(x-2/) 



-rMf, 



dr 



(47rr)^/2 



(sin 2tT)^i, 



det 



sin tJ- 



-1/2 



X exp 



1 



{x~y)- {TcoItT) ■ (x-y) 



(49) 



Second we need the vacuum polarization tensor eq.(p3|) for the two- loop effective potential eq.(p7[), whose proper time 
expression is found as follows: 



TrLn 



= ^TrLn 
2 



7^7^ (n^ - QtM^i^ - Qi^) + m"^ 



duu-'^e-""'' 



Tr 



1/A2 



where Hq has been given in eg. (|32|) , and 

Expand the final expression of eq.(|50|) with respect to the quantum field up to the second order to find 



Tr 



+ H2)\ _ 



Tr e 



-uHn 



h+h 



(50) 

(51) 
(52) 

(53) 
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with 



and 



H,{x) = Qlix) - ^<J^,,e-''Q''^{x)Ql{x)n , (55) 



U2 

2 ^0 

V? f 1 — V 1 — V 1 + V 1 + V 

— dvti- d^xd^yHi{x,y,——u)K{x,y,——u)Hi{y,x;——u)K{y,x;——u), (56) 







H,{x, y- t) = -iQ^{x) cot (x - y), - g^(x)^^,(x - y),T^ + WfQl{x)Ta + a^.Pf Q;i(a;)ra 



(57) 



where we have utihzed the \x)- representation such that 

{x\n^e-^"«\y) = (-idl + \t^,,x,T^ K{x,y;T) 



cot —] {x~ y)„I + T^^f^i'ix - y)yTz K{x, y; r) , (58) 



with the aid of eg. (|35|) , and made a change of variable from ui to v, ui = {1 + v)u/2, in the first Une of eq.(|5q). In 
terms of Iijh, the vacuum polarization tensor eq. ( ^3|) reads 

(59) 



1 *r2 poc 



If we write 



IliiA^,y)=[icosC + esmC){Ul,ix-y) + eUl,ix-y))Y' ; forz,j = l,2; (60) 



n'l(x - y) = nl^ix - y) , (Abehan Part) ; (61) 



in order to meet the expression of the gluon propagators eqs.(^) and (p^), the two- loop contribution, eq.(|27|), is 
expressed as 



V 



(D) 
2 



VT 2 



J d''xd^yIlf,{x,y)All{y,x) 

" ^ (62) 



d^p 



Again the third term designates the abelian contribution. 

The explicit forms of 11^'^^^ ~ 11^^^ are found, after performing the Fourier transformation, as follows: in 2 + 1 
dimensions, put 



0(3) = p. ( cosuvT/2~cosuT/2 \ _ 
\ uJFsinwJF/2 / ' 



(63) 

+ { cosuT l2±coBuvT I2\ ^ { sinuvT I2\ 

= -■ TFTo ; /3m- = ^ : 64 

^ \ smMj^/2 J \ smuJF/2 ' 

and utilize N^^, given in cq.(^) with obvious abbreviations such that 



= N^pNp^ ■ EE N^pNp^N^^ ; {Np),, = N^^p^ ; etc. (65) 
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to obtain 



dul dt;«i/2e-"(™'+^''') 



1 



sinhwF+/2 



P+ c -P+ / , uF+ , uvF+..^^ . F+ , uvF+.^^ 



r , , uF+ 

+ < — — — (cosh cosh 

\2sinhuF+/2^ 2 



+ F+ smh ^ ^ (iV2)^^ _ ^(cosh ^ - cosh —±){N^a-)^, 



+F+ sinh ^^(iV2/3)^, + cosh - coth ^ sinh -^6^,,p ■ {Pp) 



+ 



+ 



cosh uF+/2 + cosh / 2 
2 

cosh uF^/2 — cosh /2 



{^^.(TVp) • {Np) - {N')^.p^ +Py.{N^p). + (Af'p)^P. - 3(7Vp)^(iVp), 



+V(Ara-p) • {Na-p) - {N^Ji^^a'p) ■ (a'p) + {a' p) ^{N'^ a' p),, + {N^a-p)^,{a-p)^ - 3{Na-p)^,{Na-p),,} 



+ sinh 



uvF-L 



+ 



sinh^l - 25^,^p ■ {Na-p) + p^,{Na-p)^ + {Na-p)^,p^ - {Np)^{a-p)^ - (a"p);^(iVp)^} 



(66) 



KM = 



F~ 



dul dvu^'^e-<"''+'^'"'^ 



+ 



2(47r)3/2 7^/A2 Jo"'""' ^ sinhMF+/2 

■U-F+ _ , uvF_ 

h 3 cosh 

2 2 

cosh uF+/2 + cosh ut;F+/2 r _ i - \ \ cosh uF+ /2 — cosh uwF+/ 2 



^(cosh ^ + 3 cosh -^)N^, - ^(cosh - cosh -^){N^)^^ - F+ sinh ^(TV^a")^, 



{p^{N^a-p), 



-{N^a-p)^p. + {N^p)^{a-p)^ - {a-p)^{N^p). - 4N^.{Np) ■ (a'p) - {Np)^{Na-p). + {Na-p)^{Np)^} 

111) f 

+ sinh^|(/3p)4Arp), _ {Np)^m, + N^^p ■ {I3p)} 

-sinh^|iV,,„p2 + N^,^{a-p) ■ {a' p) + p^{Np)^ - {Np)^p^ + {a-p)^{Na-p)^ - {Na-p)^{a-p)^^ 



, ,,,, , rfu/ d^«V2e-«(-W^')_^^ 
2(47r)3/2 7,/A. Jo sinhwF+/2 



cosh uF+ /2 + cosh uvF+/2 



+ 



cosh uF+ /2 — cosh / 2 



{5ij,^p'^ -Pij,p^) +sm\i^^-^^{Np)ij,{a p)^ + {a p)^(Afp)^| 

- 5^.{Np) . {Np)-p^{N^p), - {N^p)^p. + 3{Np)^{Np),} 



+ cosh 



uF+ 



{{Pp)iJ.{Pp)t^ - Pi^i^P ■ iPp) + (a P)u^ 



Similarly in 3 + 1 dimensions, introduce 

coshuF. /2 — coshuwFi /2 , , ^. , cosh uF_ /2 — cosh 'uwF_/2 ^ , 
^ " .F^sinh.F^/2 ■ + .F_sinh.F_/2 ^''^^ ' ^^"^^ 



where 



j+ ^ F+ )iJ.v - F- . J- ^ F- )^'^ - F+ {^)nu /_ .+ ,p p 

^ijtv — pi_pi ' ^iJtv — F^ — F^ v~ ^iJ.v\^+ 

which can be regarded as projection operators, obeying 



(67) 



(68) 



(69) 



(70) 
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It + 1 

where have been defined by eq.(^). Then 



4(47r)2 



dul dt;ue-"("'+^''') ^ 



cosh — — cosh 



uvF^ 



sinh uF+ / 2 sinh / 2 

mF+ mF_ . uvF^ . uvF^ 

coth coth smh smh 

2 2 2 2 2 2 

coshuFi/2 (cosh ui^_/2 — cosh uuF_/2) , , , , , 

+^ sh.h- J_/2 -^'^'^ ■ 

cosh WW F_/2 (cosh /2 — cosh lit) F+/2) , \ /r+ ^ 

" sinh2wF+/2 ^ ^^'^^ ^^'^ 

, . uvF^ . (1 — coshui^+/2coshuuF_|_/2) (1 — coshui^_/2coshwwi^_/2) 
+ sinh — - — smh — 



F+ cosh uvF_/2 F_ cosh uF+ /2 \ 
sinhui^+/2 sinhui^_/2 / ' 

[{I+p)-{I+p)I+^-{I+p),{rp).] 



+ (+--) 



n^.(p) 



F+F^ 
4(4^2 



sinh wF_|. / 2 sinh / 2 

1 r-F- sinhui^+/2 



(7V+p)^(iV-p), 



(72) 



sinh mF+/2 sinh uF_/2 L sinhMi^_/2 

2sinhuF_/2 (coshwF+/2 — coshMwi^+/2) , 

sinh^ wF+/2 ^ ^^'^ ""^ 

uF^ , uvF^ , uvF^ 
- coth — - — smh — - — smh 



{(I+p) ■ (/+p)7V+ + {I-p)^iN+p), ~ (Ar+p)^(/-p)„} 



KM 



F+F- 



2 2 2 

coshuuF_/2(l — cosh wF+/2 cosh uuF+/2) 
sinhu,i^+/2 

:{{rp)^{N+p), - {N+p)^{rp), - iI+p)^{N+p), + iN+p)^{I+p),} 

+ (+--)] ; 

1 



(73) 



4(47r)2yi/A2 Jo 



sinh uF^ j 2 sinh uF^ j 2 



cosh 



uvF^ 



■ cosh 



uvF^ 



coth ■ 



uF^ 



■ coth ■ 



uF^ 



■ sinh ■ 



uvF^ 



■ sinh ■ 



uvF^ 



2 2 

2coshwF,/2 (coshf ,/2 - coshwwF,/2) . , _ (/+,),(/+p).} 

smh i/.h^/z I. ^ J 



{(/-p).(/-p)/+ -(/+p)^(/-p).} 



sinh^ wF+/2 

(1 - cosh wF+/2 cosh Mwi^+/2)(1 - cosh ui^_/2 cosh 'uwi^_/2) 
sinh ui^+/2 sinh uF_/2 

+ (+ 



uvFj^ . uvF^ 
— smh — - — smh — - — 



{N+p)^{N-p), 



(74) 



(These expressions are so lengthy that we ensure the correctness by checking the gauge invariance of those, that is, 
the Ward-Takahashi relation in Appendix ^.) 

Armed with these general results, in the following we consider the magnetic background only and proceed to 
calculate the gap equation in the 2 + 1 then 3 + 1 dimensions. 



III. THE GAP EQUATION IN 2 + 1 DIMENSIONS 

When the background is purely magnetic, £J = 0, in view of eq. (|39|) , F+ =^ B. A dimensionless quantity, 

is introduced in addition to x = / h? in eq. (|29|) . (Since the coupling constant e has been included to gauge fields 
the dimension of gauge fields is always one.) The one-loop contribution to the gap equation eq. (p9|) , 



9 



3/2 



reads, with the aid of eqs. (pG), (Es), and as, 



2A 9to2 



-6 



drr-i/^e-^^ coth 



(76) 



(77) 



where r has been scaled to A^t. We plot ,ff'\x) in Fig.^ It is seen that for a fixed four-fermi coupling that is, 
with respect to a (supposed) horizontal line, mass is a monotone increasing function of magnetic field strength. It 
is also noted that the critical coupling Qc, defined by f^\x = 0) in eq.(|7^) for a fixed magnetic field, goes to zero 
when B ^ 0. This phenomenon is so called "Dimensional Reduction ||l3|" and is due to the infrared divergence of the 
effective potential under the background magnetic field [E2| . 



fi (F+=B) 




0.005 



0.01, 



-3.5 



-B=0 
-B=Mg-/3 
■ B=2Mg-/3 








//' ' 








/ ' ' 


-1 


i ' 








1 1 

'! 


-2 






1 


-4 



0.5 



FIG. 1. One-loop contribution to the gap equation in 2-1- 1 dimensions. Solid line represents 5 = 0, dash-dotted line Mg'^/S, 
and dotted lines 2A4g^/3, respectively. In a:: > 0.01, recognized from the small graph, all curves become degenerate. In order 
to fix the magnitude of the background field, the (dimensionless) gluon mass y — Mg /A^ is set to be 0.01 . 

The two-loop contribution is found as. 



^2g2 .oo .1 



, ,,,,,, dui dvi ^dTu^'\-^--^y . 

2(4^)3/2A7i io io smhu6/2 



M [qi cosh 2tB + q2 sinh 2tB + q^) 
cosh tB 



-JVq4 



(78) 



where u and r have been scaled to A'^u and A^r, y, defined as 



~ A2 



(79) 



is a dimensionless gluon mass, 



M 



sinh uB / 2 



M(I-w2)/4 + r ' 



coshuS/2 — coshMwS/2 + tanh rSsinhuS/2' 



sinh uB/2 



coshwB/2 - cosh uvB/2 + tB sinh uS/2 



(80) 



and 
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JC , ^ uB , uB . , uvB. 

gi = — (cosh — V coth — smh ) 



92 
94 



^, ,, , uB , uuS cosh — cosh (1 — coshuS/2 coshuuS/2) 

-B M cosh — cosh ' i ^ 



smh^ uB/2 



/C , . , uB 
2(«inh — 



. , uwS, „ / . coshiiuS/2 , , uB 
V smh — — ) + B [ M —r-, — TTTTT (cosh — 
2 ^ V sinhw6/2 ^ 2 



cosh ) — cosh 



sinh uB/2 
uvB^ 



IC n uvB 

— (V cosh— — 
4 ^ 2 

IC . , uvB 

— (cosh — — 
4 ^ 2 



, , uvB. ^ , cosh uS/2 — cosh uwS/2 1 / , uvB 
V coth — smh — — ) + BM — ^— \ — I cosh ■ 



sinh^ uB/2 



uB 



2smh'uS/2 



uB . uvB . 

V coth — smh 1 

2 2 ^ 



JV / cosh uB /2 — cosh uvB/2 cosh uwS /2 — w coth W;B/2 sinh itwS / 2 



sinh^ Me/2 



(81) 
(82) 
(83) 

(84) 



qi,q2, and terms come from the first two terms in cq.(|62|), that is, from the nonabeUan contribution. Meanwhile 
54 represents the abehan contribution. 



0.003 



0.005 



-0.25 




-0.35 



FIG. 2. Two-loop contribution to the gap equation in 2 + 1 dimensions. As in Fig.l, the same line pattern is used. From 
the smaller graph, we see that all curves become degenerate but remain negative in the whole region 1 > a: > 0. The graph is 
drawn by putting e^/(47rA) = 0.01 ,y = 0.01. 



We plot the two-loop part in Fig.y. (Our choice of the gauge coupling e^/(47rA) — 0.01 does guarantee the 
approximation; since / ~ 0.1 by comparing the vertical scale between Fig.|^ and Fig.^.) From the small graph, 
we can convince that gluons enhance x^B as is expected; because all curves remain negative for a whole region, 
< a; < 1. Note that there is a crossover around x ~ 0.001: in the region larger than the crossover, a;, on some 
horizontal line (a line with a fixed four-fermi coupling), is a monotone decreasing function of the magnitude of the 
background field. In the region smaller than that, x is, however, a increasing function of it, similar to the one-loop 
case. To see the situation more carefully, we plot the abelian contribution to the gap equation, that is, q^ term in 
eq.(|7S 
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0.001 



0.003 



0.005 




B=0 

B=Mg^/3 
B=2Mg2/3 



-0.15 



FIG. 3. Abelian part of the two-loop contribution to the gap equation in 2 + 1 dimensions. The graph is drawn by putting 



e7(47rA) = 0.01 ,y = 0.01. 

From Fig.^, x, on some horizontal line, is a monotone decreasing function of the background field everywhere for a 
fixed g. Therefore, the increasing tendency of Fig.H in a; < 0.001 comes from the nonabelian parts in eq.dTq). 



IV. THE GAP EQUATION IN 3 + 1 DIMENSIONS 

In 3 + 1 dimensions, when E = 

F+=^ B , F-=^0 . (85) 
Again employing the dimensionless quantity B = B / h?, we have the gap equation of one-loop contribution, 

4(7r)2 



with 



^ '~ 2A2 9m2 



tB 



-B I drr-^e-^^coth 
li 2 



(86) 



(87) 



which is plotted in Fig.^. All curves become degenerate again where x is large. For a fixed four-fermi coupling g, 
mass is a monotone increasing function of magnetic field strength. Moreover the critical coupling goes to zero even 
under infinitesimal magnetic fields [13 12 1. 



12 



f 1 (F+=B, F.=0) 



0.001 



0.003 



0.005, 



1.95 



-B=0 
-B=Mg-/3 
B=2Mg-/3 




FIG. 4. One- loop contribution to the gap equation in 3 + 1 dimensions. Solid line, dash-dotted line, and dotted line designate 
5 = 0, Mg /3, and 2Mg /3, respectively. The smaller graph is shown the whole structure, 1 > a; > 0. We again set M^/A^ = 0.01. 



The two-loop contribution is given by 



^2g2 .oo ,1 



2(47r)2yi Jq J„ sinhu6/2 _ 



where K.,M, and N are the same as eq.(pO[) and 
pi = /C(l - w^) cosh 



M {pi cosh 2tB + P2 sinh 2tB + ps) 



cosh tB 



-NpA 



9n , uB „ . , f , uvB 
^ "-""^ + BM cosh 



uB . uvB cosh U;B/2 — cosh 

V coth smh ^ 

2 2 -"^2 



sinh uS/2 



B 



2 / , 

- — cosh 

u\ 2 2sinhii6/2 



. , ui3 ^ . /cosh wS/2 cosh iiuS/2 — 1 . , 2sinhwS/2 
P2 = JC{l-%r)s\nk—+BM[ —. — t;sinh- ' 



P3 



uvB 



K. ( cosh 
V 2 



u coth — smh 

2 2 



sinh uS/2 
1 



cosh 



uvB 



2BM 



coshwS/2 — coshuwS/2 
sinh^uS/2 
2 / uB 
\^ ~2 2sinhw;B/2 



Pi 



^ ( cosh „ 
2 V 2 



uvB 



uB . uvB 

V coth smh 

2 2 



1-^2 

— - — cosh — — 
2 2 



— (cosh 
2 V 



itwS uuS coshwS/2 — coshuwi3/2 

V coth — smh 1 ^ 

2 2 2 sinh2uS/2 



(89) 
(90) 

(91) 

(92) 



The graph is shown in Fig.^ (The choice of the gauge couphng e'^/Air = 0.01 again guarantees our approximation; 



f(4) / .(4) 

Yamawaki [g2[.) All curves, shown in the smaller graph, remain negative and become degenerate for x > 0.001 



since ff>/n^> - 0.05. It should be noted that our result of S = is consistent with that of Kondo, Shuto, and 
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J, (4) 

f 2 (F+=B. F.=0) 



0.0005 



-0.099 




0.001, 



-0.0995 



FIG. 5. Two-loop contribution to the gap equation in 3-1- 1 dimensions. In the smaller graph the whole structure, < a; < 1, 
is shown. The same line pattern is used for different curves by putting e^/47r = 0.01 and / = 0.01. 



Therefore, gluons enhance xSB everywhere even in no background _B = 0, which fits our expectation. Contrary to 
the case in 2 -I- 1 dimensions, mass is a monotone increasing function of magnetic field strength everywhere for a fixed 
four-fermi coupling g. 



V. DISCUSSION 



In this paper we discuss the effect of dynamical SU{2) gluons to the gap equation of NJL model under the influence 
of the constant background magnetic field. The two-loop calculations make expressions considerably complicated but 
correctness of the results is guaranteed by checking the Ward-Takahashi relation in Appendix In 3 + 1 dimensions, 
as is seen from Fig.||and Fig.^, gluons play the same role as fermion in the one-loop, that is, enhance x^B. Moreover 
dependence of gluons on the background field is also same as fermion in the one-loop: dynamical mass grows larger as 
the background magnetic field becomes stronger. The result is consistent with our expectation but different from 
where RG with the one-loop calculation was employed. In 2 -I- 1 dimensions, the situation is unchanged that gluons 
enhance x^B even under influence of the background field, contrary to the work of [T7| where a different gauge was 
employed. Dependence of gluons on the background magnetic field, however, is not so simple as in 3 -f 1 dimensions: 
as is seen from FigJ^, when dynamical mass is tiny the background field increases it, but in a well-broken region, that 
is, in a region where dynamical mass is large, the background field resists a mass to grow. Difference between 2 + 1 
and 3 + 1 dimensions is due to that of the u dependence in eqs.(|7^) and (p8|): by making the scale transformation to 
u and T such that 



u 
B 



T 

B 



quantities, IC,M,Af (|80|), scale 



so that Qi and pi 



Therefore, 



K^BJC: M^M: Af ^ N 



1^4), 



) and (89) ~ (pa) transform 



Bqi ; ^ Bpi 



(93) 



(94) 



(95) 





VBf['^ 




- Bfi'^ 


f(4) 
J2 





(96) 
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which shows a monotone character with respect to B, convincing us of the results of Fig.|i| , Fig.^ , and Fig||, 
quahtatively. (Because there is still B dependence in f[^\ /{^'' , and f2'^\) But 



(3) 



(3) 



(97) 



which implies that dependence on B is due only to the detailed structure of the integrand of the expression (|78|), that 
is, we cannot extract a simple monotone behavior in this case. 

The second point we wish to discuss is on instability of giuon functional determinant: we have avoided this by 
introducing gluon mass Mg which is assumed always bigger than the magnitude of background magnetic field B; 
Mg > B. Physically, it is interpreted that the energy of nonabelian particles could become lower and lower as the 
background magnetic field grows larger and larger. There is no lower limit in the system. The situation is exactly 
the same in the constant electric field case, where the vacuum becomes unstable due to successive pair productions. 
We have treated this pathological instability by considering only external electric field whose magnitude E is less 
than that of dynamical mass squared; vn? > E The point is that the setup itself- "field theories under constant 
background field"- is pathological. The system is not closed; energy is continuously supplied from outer environment. 
However, even in these pathological environments, we could still think about those background effects to the systems 
provided that their magnitude is so small. 

The final point to discuss is beyond the tree approximation of the auxiliary fields a and tt: in most cases of NJL- 
study, fermions are assumed to have N components with N being supposed infinite finally. However in the actual 
situation N is finite so that 0{1/N) corrections should be taken into account. A study in a simpler model |^ says 
that the approximation becomes more and more accurate if wc incorporate higher order terms. Thus going beyond 
the 1-loop of the auxiliary fields is captivating and the work in this direction is under progress. 



APPENDIX A: CALCULATION OF KERNEL BY THE PROPER TIME 



In this appendix we derive the expression of the kernel (p4|). 



(Al) 



Because of the covariantly constant condition (||), the matrix element {x\e ^'^f [y) can be calculated exactly the same 
way as the abelian case ; 



(4^t)^/ 



det 



sin tT/2 



exp 



1, , tT\ , 



(A2) 



with 



C — ^yu 



(A3) 



The remaining task is therefore the calculation of exp 
given as 



2 ^ iiy^'i 



in 2 + 1 dimensions, the gamma matrices are 



7m 



CJu 



-(Tu 



Cut/ _ 1 



<Jn 



(A4) 



where J^'s satisfy, 



(A5) 



In terms of Ju's, 



-cr^yJFpy — E2J1 — E1J2 + BJj, ; B — T\2 ; E — (.F13, JF23) 



(A6) 



From eqs. (|A5|) and (A6), we obtain 
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Meanwhile, 



Therefore in 2 + 1 dimensions, 



exp 



tF+ tF+ 
= cosh —^1 + (Tf^„Nf_,^ smh — 



In 3 + 1 dimensions, first write 



\\B + E\ + IB - E\} 
_ {\B + E\-\B-E\} 



and introduce the antisymmetric tensors, N^^ = — -/Vj^, such that 



N- = 



p2 _ p2 ' "'pi' 



Fl-Fl 



which satisfy 



where the second relation can be verified by using eqs.(AlO) and 

2 



With the aid of Nj^^, T^j^v is expressed as 



givmg 



By noting 



and eq.(|AT^), we find 



J'^, = F+N+, + F-N^^ 



Therefore 



Also by noting 









exp 




= exp 



F^ 



exp 



F^ 



(A7) 



(A8) 



(A9) 



(AlO) 



(All) 

(A12) 
(A13) 

(AM) 

(A15) 

(A16) 

(A17) 
(A18) 

(A19) 



(A20) 
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and eq.( A13| ), 



Hence 



yielding 



Finally utilizing 



we obtain 



exp 



± \2 



exp 



= cosh + f^t^uN^^^ sinh 



C^M^^^^crpAA^pA = -475 



2 ^^U'^ ui^^Z 



cosh — — cosh 75 sinh — — sinh ) I 

2 2 2 2 



/ TF4. tF^ tF^ tF- 

,ylN+^ sinh cosh + A^^^ cosh sinh ) T3 



flu y ■ 2 

= Xo(r)I + if3(T)T3 : 

Ti^4. rF_ tF+ tF^ 

Ko{t) EE cosh cosh — ^ 75 sinh sinh — ^ , 

Kz (r) EEE CTp^ ( sinh — - cosh — ^ + N^^ cosh — - sinh — 



APPENDIX B: GLUON PROPAGATOR IN TERMS OF THE PROPER TIME 

In this appendix we show the proper time representation of the gluon propagator 

where we have introduced a slightly different notation from eg. (pi]), 






The proper time representation is obtained as usual; 





(x|e--n^|y) 


'0 





rjr[ad(T3)] 



cb 



where 

EE {nff ■ EE - ^^(i)[ad(T3)]"^ , 
and we have introduced the gluon mass Mg to avoid the tachyonic singularity. For a = 6 = 3, it reads 

which is the free propagator. Therefore we obtain eq.(p7|): 

1 



A' (a;-y) = 5,,. / dr 



(47rT)^/2 



exp 



4t 
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In a, 6 = 1, 2 use i, j = 1, 2 and utilize the result in Appendix ^ (|A2|) and [ad(T3)]*-'' = -ie'^ to find 



\y) 



y 1 



r(^''^cosC + e*^sinC) 



det 



sin tJ- 



exp 



(47rT)^/2 

where C is given in eq.( [A3| ). Finally by noting that 

/ ~2..^[ad(T3)]y^ =<5'^(cos2T^)^,-e'^"(sin2r^)^. , 



(a; - 2/)^ (J^cotrJ^)^^ (a; - y). 



, (B7) 



(B8) 



the relations (|46D ~ are obtained 



APPENDIX C: PROOF THAT OUR CLASSICAL SOLUTION SATISFIES THE COVARIANTLY 

CONSTANT CONDITION 

In this appendix we show that the right-hand side of eq.(p^ vanishes: 



= 0, 



tr (7p(5p - i^^Ts) + to) ( - Tq) = tr 5*^(0;, x) ( - i^yTa) 
where Sa{x, y) is the fermion propagator under the background fields; 

(7m (^M ^ i-^^Ta) + ■m)SA{x,y) ^ 5° {x - y) . 
Sa{x, y) can be expressed, by using proper time method, as 

SA{x,y)^ {x\{i-i^Ii^+m) ^\y) = (a:;|(-i7^n^ + m)(n^ - iapA^pATs + to^)"^ |y) 
dre-™' {x\ (-^7M^M + m)e--«° \y) 



dre"™' ^i7p^^cot^^ {x - y)J. - ]p^T^,y{x - y)i,Tj, + rri\ K{x,y-,T) , 



(CI) 



(C2) 



(C3) 



where K{x, y; r) is the kernel of eq.( Al). Therefore 



tr 



Sa{x,x){ - i-iuTa) 



'0 



det 



sinrJ^/2 
tT/2 



fJ.V 



-1/2 



tr [ei'^'-^^'-^^'j^Ta] = ; (C4) 



since the trace for the gamma matrices vanishes because the total number of those is odd. 



APPENDIX D: THE WARD-TAKAHASHI RELATION OF VACUUM POLARIZATION 

In this appendix it is shown that the vacuum polarization function satisfies the Ward-Takahashi relation, 



where 



with being a background field 



iV;r'llll{x,y)^0, 



(Dl) 



The Ward-Takahashi relation (Dl) is separated into 



(D2) 
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(D3) 



and 



(^5^1 _ e'i^A^,{x))Ilf,{x, y) = ; for k^l,2 



(D4) 



The first relation (D3) can easily be checked by noting eq.(|68|) and eg. (|74|) so that the second relation (D4) must be 
examined. In view of the fact that TV^j^{x^y) can be written as 



^%{x,y) = [(cosC + esinC) {Wl^{x - y) + eW%[x ^ y))\ 



the matrix relation reduces to 



dinl,{x -y) + A^{x- y)n'^,{x - y) = 

dn^lA^ -y)-A^{x- y)ni (a: -y) = {) 



(D5) 

(D6) 
(D7) 



Utilizing the series expansion with respect to the background gauge field, we show, up to 0{F), these relations indeed 
hold: first note that 11^^ and 11^^, are polynomials of even and add powers of F respectively. Thus in 0(1) eg. (p^) 
reads 



F=0 



= , 



which is fulfilled; since from eq.(|6^) and eg.(|7^), 11^^ has been given by 

1 



Next in 0{F), 



^=0 (4^)^/2 



1 - t;2 
u 



du I dv 



D/2-1 "^^P 



2 1 - 2 

u \ ni H p 



0(F) 



F=0 



which becomes in the momentum space, 



F=0 



= , 



where (p) 



0(F) 



can be found from the expression (|67|) in 2 + 1 dimensions. 



KM 



0(F) (47r)3/2 



du / dvu 



exp 



u [ m + 



2 , 1-^^ 2 



and from eg. ( |7^ ) in 3 + 1 dimensions 

KM) ^ 



, „ , du I dv exp 
o{F) (47r)2 y^/A^ Jo 



, 2 1 - 2 

TO H p 



<! -N, 



({rp)^{N+p), ~ (iV+p)^(/-p), - {I+p)^(N+p), + (7V+p)^(/+p),) 

2 ^ 

+y • + {I-p)^{N+p\ - {N+p)^{I-p\) [ + (+---) 



(D8) 



(D9) 



(DIO) 



(Dll) 



(D12) 



(D13) 



respectively. With the use of eg.(|71|), the left-hand side of eg. (Dll) is shown to vanish. 
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L.H.S 



1 ^—um' 

du I dv 



dv 



v{l — v'^) exp 



-P 



= . 



(D14) 



Therefore we can convince ourselves that the Ward-Takahashi relation is satisfied in each order of the background 
field. 
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